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Abstract 

A model for the dark halos of galaxy clusters, based on the Weyl 
geometric scalar tensor theory of gravity (WST) with a MOND-like 
approximation, is proposed. It is uniquely determined by the baryonic 
mass distribution of hot gas and stars. A first heuristic check against 
empirical data for 19 clusters (2 of which are outliers), taken from 
the literature, shows encouraging results. Modulo a caveat resulting 
from different background theories (Einstein gravity plus ACDM ver¬ 
sus WST), the total mass for 14 of the outlier reduced ensemble of 17 
clusters seems to be predicted correctly (in the sense of overlapping 1 a 
error intervals). 


Contents 

Introduction 

1 Theoretical framework |4] 

1.1 Weyl geometric scalar tensor theory of gravity (WST). |4] 

1.2 The weak field static approximation . IH] 

1.3 WST gravity with cubic kinematic Lagrangian (WST-3L). [10] 

1.4 A WST approach with MOND-like phenomenology . |T3| 

1.5 Comparison with usual MOND theories. Uni 

1.6 Short resume. [Ill 

2 Halo model for clusters of galaxies 1201 

2.1 Cluster models for baryonic mass (hot gas and stars) . [501 

2.2 Two contributions to the scalar field energy in clusters of galaxies . . |2T] 

2.3 Scalar field and phantom halos Psf i, Pphi in th® cluster-barycentric 

MOND approximation. 

2.4 Scalar field halos of galaxies in their respective galacto-centric MOND 

approximations. [^ 

2.5 A three-component halo model for clusters of galaxies. |23| 


’University of Wuppertal, Faculty of Math. & Nat. Sciences and Interdisciplinary 
Centre for Hist, and Phil, of Science; scholz@math.uni-wuppertal.de 


1 













3 A first comparison with empirical data 1251 

3.1 Empirical determination of mass data for galaxy clusters . 

3.2 Theory dependence of mass data for galaxy clusters . UHl 

3.3 Empirical data for 19 clusters . UHl 

3.4 Adapting the WST halo model to the data . HH 

3.5 The WST halo model with the Coma cluster as test case . . .... ED 

3.6 Halos and total mass for 17(+2) clusters of galaxies . EH 

3.7 Comparison with TeVeS and NEW halos . EHl 

3.8 A side-glance at the bullet cluster . EOl 

4 Discussion I4T] 

5 Appendix 1431 

5.1 Energy component of the scalar field. EH 

5.2 Remarks on the numerical implementation. EH 

Bibliography 1451 


Introduction 


This is a CORRECTED, and in some lesser aspects improved, POSTPRINT ver¬ 
sion of the paper of the same title, published in JOURNAL OF GRAVITY (JG) 
volume 2016, article ID 9706704- The most important corrections refer to 
equs. (28) and (47) of the publication in JG. An error in JG (47) lead to 
a wrong relation between the acceleration due to the scale connection (a^p) 
and the acceleration arising from the scalar field energy density (ogf). The 
corrected form of JG (47) appears here as equ. (46) and has consequences 
for the model (equs. (50)ff.). It makes a new run of the data evaluation 


necessary, now based on the corrected dynamical equations. The new results 
are given in the updated tables 2, 3 and figs. 2 to 7. The overall picture of 
the empirical test does not change, although now three rather than two galaxy 
clusters agree with the model only in the 2a range. 


In this paper the gravitational dynamics of galaxy clusters is investigated 
from the point of view of Weyl geometric scalar tensor theory of gravity 
(WST) with a non-quadratic kinematic Lagrange term for the scalar field 
(3L) similar to the first relativistic MOND theory rAQUAL (“relativistic a- 
quadratic Lagrangian”) [23]. To make the paper as self-contained as possible, 
it starts with an outline of WST-3L (section 1). WST-3L has two (inho¬ 
mogeneous) centrally symmetric static weak field approximations: (i) the 
Schwarzschild-de Sitter solution with its Newtonian approximation, which is 
valid if the scalar field and the WST-typical scale connection plays a negli¬ 
gible role; (ii) a MOND-like approximation which is appropriate under the 
constraints that the scale connection cannot be ignored but is still small 
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enough to allow for a Newtonian weak field limit of the (generalized) Ein¬ 
stein equation. The acceleration in the MOND approximation consists of 
a Newton term and an additional acceleration of which three quarters are 
due to the energy density of the scalar field and one quarter to the scale 
connection typical for Weyl geometric gravity. 

In centrally symmetric constellations the scalar field energy forms a halo 
about the baryonic mass concentrations. Besides the acceleration derived 
from the (Riemannian) Levi-Civita connection induced by the baryonic mat¬ 
ter and the scalar field energy an additional acceleration component due to 
the Weyl geometric scale connection arises in the present approach. If the 
latter is expressed by a fictitious mass in Newtonian terms, a phantom halo 
can be ascribed to it. It indicates the amount of mass one has to assume 
in the framework of Newton dynamics to produce the the same amount of 
additional acceleration. The scalar field halo consists of true energy derived 
from the energy-momentum tensor; it is independent of the reference sys¬ 
tem, as long as one restricts the consideration to reference systems with low 
(non-relativistic) relative velocities. The phantom halo, on the other hand, 
is a symbolical construct and valid only in the chosen reference system (and 
scale gauge). For galaxy clusters we find two components of the scalar field 
halo, one deriving from the total baryonic mass in the MOND approximation 
of the barycentric rest system of the cluster (component 1), and one arising 
from the superposition of all the scalar field halos forming around each single 
galaxy in the MOND approximation of the latter’s rest system (component 
2). Because velocities of the galaxies with regard to the cluster barycenter 
are small (non-relativistic) and also the energy densities are small, the two 
components can be superimposed additively (linear approximation). With 
regard to the barycentric rest system of the cluster a three component halo 
for clusters of galaxies arises, two components being due to the scalar field 
energy and one purely phantom (section 2). 

The two component scalar field halo is a distinctive feature of the Weyl 
geometric scalar tensor approach; it is neither present in the non-relativistic 
MOND approaches nor in rAQUAL. One may pose the question whether it 
suffices for explaining the deviation of the cluster dynamics from the New¬ 
tonian expectation without additional dark matter. If we call the totality of 
the three components the transparent halo of the cluster (section 2.5), the 
question is whether the (theoretically derived) transparent halo can explain 
the dark halo of galaxy clusters, observationally determined in the framework 
of Einstein gravity and ACDM. 

Section 3 contains a first test of the model by confronting it with em¬ 
pirical data on mass distribution available in the astronomical literature. A 
full-fledged test would presuppose an evaluation of raw observational data in 
the framework of the present approach and is beyond the scope of this study. 
Here we use recently published data on the total mass (dark plus baryonic), 
hot gas, and the star matter for 19 galaxy clusters, which have been deter- 
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mined from different observational data sources (and are thus more precise 
than earlier ones) |27) . [25], |28] . 2 of the 19 clusters show a surprisingly 
large relation of total mass to gas mass. They are separated as outliers from 
the rest of the ensemble already by the authors of the study; so do we. 17 
non-outlying clusters remain as our core reference ensemble. 

In the mentioned studies total mass, gas mass and star mass are deter¬ 
mined on the background of Einstein gravity plus ACDM. That raises the 
problem of compatibility with the WST framework. It is discussed in 3.1, 
3.2 and leads to a certain caveat with regard to the empirical values for the 
total mass (M 20 O) ^hoo) and the reference distances r 200 ) ^500 to the cluster 
centers. But it does not seem to obstruct the possibility for a hrst empirical 
check of our model (section |3.2| ). More rehned studies are welcome. They 
have to use the WST framework for evaluating the observational raw data or, 
at least, to analyze the transfer problem of mass data from one framework 
to the other in more detail. 

For 14 of the 17 main reference clusters the empirical and the theoretical 
values for the total mass agree in the sense of overlapping 1 a error intervals. 
For the remaining three we hnd overlap in the 2 a range. The two outliers of 
the original study do not lead to overlapping intervals even in the 4 a range 
(section 3.6). In the present approach the dynamics of the Coma cluster is 
explained without assuming a component of particle dark matter. It is being 
discussed in more detail than the other clusters in section 13.51 

A short comparison with the halos of R. Sanders’ /r 2 -MOND model with 
an additional neutrino core |20| . and with the NFW halo [13] is given for 
Coma (section |3.7[ ). The paper is rounded off by a short remark on the bullet 
cluster (section |3.8[) and a hnal discussion (section Q . 


1 Theoretical framework 

1.1 Weyl geometric scalar tensor theory of gravity (WST) 

Among the family of scalar tensor theories of gravity the best known ones, 
and closest to Einstein gravity, are those with a Langrangian containing a 
modihed Hilbert term coupled to a scalar held (f). Their Lagrangian has the 
general form 

L = + (1) 

= Ly/\^\, \g\ = \detg\. 

Here g is an abbreviation for a 4-dimensional pseudo-Riemannian metric 

9 = (ff/ii/) of signature (-h -|-+). (p is a real valued scalar held on space- 

time, its kinetic term, ^ a constant allowing for bridging the hierarchy 
between the Planck scale and cosmologically small quantities (like vacuum 
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energy density), and the dots indicate matter and interaction terms. Under 
conformal rescaling of the metric, 

diiu g'fiu — 9fj-v ^ positive real valued function), (2) 

the scalar field changes with weight —1, i.e. (f) ^ (f)' = (/>. So far 

this is similar to the well known Jordan-Brans-Dicke (JBD) scalar-tensor 
theory of gravity M: 0, E], |in| . But here we work in a scalar-tensor 
theory in the framework of Weyl’s generalization of Riemannian geometry 

[aiisiiiTlEaEa]. 

Crucial for the Weyl geometric scalar tensor approach (WST) is that the 
scalar curvature R and all dynamical terms involving covariant derivatives 
are expressed in Weyl geometric scale covariant form0 The Lagrangian den¬ 
sity C is invariant under conformal rescaling for any value of the coefficient 
of the modified Hilbert term. For the matter and interaction terms of the 
standard model of elementary particles, scale invariance is naturally ensured 
by the coupling to the Higgs field which has the same rescaling behaviour as 
the gravitational scalar field. Although there is no complete identity, there 
is a close relationship between scale and conformal invariance in quantum 
field theory |15| . For classical matter we expect that a better understanding 
of the quantum to classical transition, e.g. by the decoherence approach, 
allows to consider scale invariant Lagrangian densities also. For the time 
being we introduce the scale invariance of matter terms in the Lagrangian 
as a postulate. 

In our context an important consequence is the scale covariance of the 
Hilbert energy momentum tensor 

^ _ 2 6Cm 

which is of weight w{T^jf) = — 2. That is consistent with dimensional con¬ 
siderations on a phenomenological level. It has been shown that the matter 
Lagrangian of quantum matter (Dirac field, Klein Gordon field) is consistent 
with test particle motion along geodesics (autoparallels) 7 (r) of the affine 
connection, if the underlying Weyl geometry is integrable (see below, equ. 
([^) [2]. For classical matter we assume the same. It is to be expected that 
it can be proven similar to Einstein gravity. 

We do not want to heap up too many technical details; more can be found 
in the literature given above. But we have to mention that a Weylian metric 
can be given by an equivalence class of pairs (g, ip) consisting of a pseudo- 
Riemannian metric g = g^udx^dx’^, the Riemannian component of the Weyl 

^Fields X are scale covariant if they transform under rescaling by X i—>■ A = Q,™X, 
with ui G Q (in most cases even in X); w is called the weight of X. Covariant derivatives 
of scale covariant fields are defined such that the result of covariant derivation is 

again scale covariant of the same weight w as X. 


(3) 
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metric, and a differentiable one-form ip = ip^dx^, the scale connection (or, in 
Weyl’s original terminology, the “length connection”), is often called the 
Weyl covector or even Weyl the vector held. In fact, p denotes a connection 
with values in the Lie algebra of the scale group , •) and can locally be 
represented by a differentiable 1-form. The equivalence is given by rescaling 
the Riemannian component of the Weylian metric according to ([^, while 
p has the peculiar gauge transformation behaviour of a connection, rather 
than that of an ordinary vector (or covector) held in a representation space 
of the scale group: 

^1, = - -jf ( 4 ) 

or, shorter, p' = p — dlogfl. 

A Weylian metric has a uniquely determined compatible affine connection 
r; in physical terms it characterizes the inertio-gravitational guiding held. 
It can be additively composed by the well known Levi-Civita connection 
^ of the Riemannian component g of any gauge {g, p) and an additional 
expression in the scale connection, in short 


r = ffT -hipT with = 5>i;px + 5^pu - gu\P^- 


(5) 


Covariant derivatives D in the Lagrangians (Q, (below) ( |28[ ), and conse¬ 
quently in the expression for the energy-momentum tensor of the scalar held 
(35 36) below, denote those of Weyl geometry. For a covariant held of 
weight w the derivativation V with regard to F of (|^ H is supplemented by 
a term due to the scaling weight w oi X: 


D^X^ = + w p^x^ = d^x^ + T^y^X^ + w p^xr (6) 

It turns out that for the metric the full covariant derivative is zero: 


D\gfj,i/ — 0 i > Xxg^i^ “^.pxdtiu — 0 (7) 


This is the Weyl geometric compatibility condition between metric and affine 
connection (sometimes called “semi-metricity”). 

In the low energy regime there are physical reasons to constrain the scale 
connection to the integrable case with a closed differentiable form dp = 0, 
i.e. d^pi, = dyp^. Then is a gradient (at least locally) and may be given 
by 

. ( 8 ) 

This constraint is part of the dehning properties of WST. Then it is possible 
to “integrate the scale connection away” |2]. Having done so, the Weylian 
metric, given as {g,p), is characterized by its Riemannian component g^i, 
only (and a vanishing scale connection = 0). By obvious reasons we call 
this the Riemann gauge. This is the analogue of the choice of Jordan frame 
in JBD theory. 
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In any case, the choice of a representative also fixes both together 
define a scale gauge of the Weylian metric. Conformal rescaling of the metric 
is accompanied by the gauge transformation of the scale connection Q. 
From a mathematical point of view all the scale gauges are on an equal 
footing, and the physical content of a WST model can be extracted, in 
principle, from any scale gauge. One only needs to form a proportion with 
the appropriate power of the scalar field.From the physical point of view 
there are, however, two particularly outstanding scale gauges. Of special 
importance besides Riemann gauge is the gauge in which the scalar field is 
scaled to a constant ipo {scalar field gauge). For the particular choice of the 
constant value such that 


{Cfi,)^ = {87rG)-^ = E^^i, ( 9 ) 

with the Newton gravitational constant G, this gauge is called Einstein gauge 
{Epi the reduced Planck energy). It is the analogue of Einstein frame in 
JBD theory. In this gauge the metrical quantities (scalar, vector or tensor 
components) of physical fields are directly expressed by the corresponding 
field or field component of the mathematical model (without the necessity 
of forming proportions). 

In Riemann gauge F reduces to ^F by definition. Thus, in this gauge, the 
guiding field is given by the ordinary expression of the Levi-Civita connec¬ 
tion. On the other hand, in Einstein gauge the measuring behaviour of clocks 
are most immediately represented by the metric field; and also other physi¬ 
cal observables are most directly expressed by the field values in this scale. 
Then the expression of the gravitational field and with it the expressions for 
accelerations contain contributions from the Weylian scale connection. Thus 
a specific dynamical difference to Einstein gravity and Riemannian geometry 
(as well as to JBD theory) arises even in the case of WST with its integrable 
Weyl geometry. 

Writing the scalar field in Riemann gauge {g, 0) in exponential form, 
(j) = e'^, turns its exponent 

w := ln0 (10) 

into a scale invariant expression for the scalar field. In the following we shall 
omit the tilde sign to simplify notation. The scale connection ip = in scalar 
field gauge is then 

(f = —duj , (11) 

because D = ^ = is the rescaling function from Riemann gauge to scalar 
field gauge. For more details see [Ill6l[l6llI7l|23]. 

For the sake of consistency under rescaling we consider scale covariant 
geodesics ^{t) with scale gauge dependent parametrizations of the geodesic 
curves of weight w{fi) = — 1: 

— ipp^u^u^ = 0 , = 7^ (12) 
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Here the affine connection contains a (/^-dependent term in addition to the 
well-known Levi-Civita connection derived from (eqn . The last 
term on the l.h.s. of ( 12 ) takes care for the scale dependent parametrization 


(compare (|^). In this way we work with a projective family of paths. 

For any gange of the Weylian metric and the scalar field, {g,ip,(f)), any 
timelike geodesic has thns a generalized proper time parametrization 7 (r) 
with g^yU^vy = — 1, where = 7 ^. Inverting the coordinate time fnnction 
t(r) along the geodesic by r(t) we have, in abbreviated notation, r't' = 1 
and thns: 


dr d dr dx® / dr 
dt'^ dt dr dt dr I dt 




dr^ 


-(^) (13) 


dH dx® 
dr^ dr 


With (eqn. (|12[)) and indices i,j, k = 1, 2,3, ... = 0,1, 2, 3 this leads to 


d^x® 

dt'^ 


fdrV f . dx^dx®® dx® dx^ 

(s-j (-Tiirir+ 


dr \ dx® 


^ dt j dr I dr dr 


dr dr 
dx^ dx®® dt dx^ 


+ ‘Pfj- 


■^0 


dr dr 
dx^ dx®® dx® 


(14) 


j dx^ dx®® dx® dx^ , 

^®® dt dt dt dt ^®® dt dt dt 




dx^ dx'' 
dt dt 


Happily, the length connection terms coming from the scale covariance mod¬ 


ification of the geodesic eqnation (12) cancel. This is an expression of the 


fact that only the trace of the geodesic enters into the coordinate time 


parametrization of the dynamical eqnation (15). The eqnation of motion 


for mass points in Weyl geometric gravity, parametrized in coordinate time, 
becomes: 






= -r; 


dx® 


n I pO 

oo + ioo^ 

dx' dx^ «X” ax-' ax' 


2r — 

dt 

dx' dx^ dx^ 


dx^ dx^ 


dt 


(15) 


In the resnit the dynamics of mass points in Weyl geometric gravity is gov¬ 
erned by the gniding field (the affine connection), like in the semi-Riemannian 
case 


eqn. (9.1.2)]. Note, however, that in (15) the length connection en¬ 
ters into the ajjine connection and inflnences the dynamics becanse of ([^. 

The geodesic eqnation thns contains terms in the scale connection In 
the low velocity, weak field regime the eqnation of motion rednces to the form 


well known from Einstein gravity = —Tqq. Here the Tgg = gTgg -t-i^cTgg 
(j = 1, 2, 3) are the coefficients of the Weyl geometric affine connection with 
<^;^Q given by (j^. This is the crncial modifying term for gravity in the Weyl 
geometric approach (low velocity case). 
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1.2 The weak field static approximation 


We want to understand the scale connection for the motion of point particles. 
The free fall of test particles in Weyl geometric gravity follows scale covariant 
geodesics. It is governed by a differential equation formally identical to the 
one in Einstein gravity (15). Here we look at the weak field static case 


for low velocities in order to study the dynamics of stars in galaxies and 
galaxies in clusters. For studies of the gas dynamics and its modification in 


our framework the velocity dependent terms of (15) have to be taken into 
account. This is not being done here. 

Analogous to Einstein gravity, the coordinate acceleration a for a low 
velocity motion x{t) in proper time parametrization is given by 




-W . 

^ oo • 


According to ([^ the total acceleration decomposes into 

= -g^io 

(j = 1, 2, 3 indices of the spacelike coordinates), where 


(16) 


(17) 


~ g^oo 


(18) 


is the Riemannian component of the acceleration known from Einstein grav¬ 
ity. Clearly 

< = (19) 

represents an additional acceleration due to the Weylian scale connection. 
For a diagonal Riemannian metric g = diag {goo, ■ ■ ■, 933 ) the general expres¬ 
sion ^ simplifies to — 1^00 = —goo9^^- General considerations on observable 
quantities and consistency with Einstein gravity show that, in order to con¬ 
front it with empirically measurable quantities, we have to take its expression 
in Einstein gauge if we want to avoid additional rescaling calculations |221 
sec. 4.6]. 

For a (diagonalized) weak field approximation in Einstein gauge. 


gfiu — T ^ 1) (^9) 


with r] = Csig diag( —1, -)-l, +1, +1), the Riemannian component of the accel¬ 
eration is the same as in Einstein gravity. Its leading term (neglecting 2-nd 
order terms in h) is. 


J — ri ~ 

“r ~ 9-^ 00 

-rj^^djhoo (no summation over j) . 

(21) 

In the limit, '■= —\hoo 

behaves like a Newtonian potential 



Or Ri — V<l>Ar, 

(22) 
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where V is understood to operate in the 3 spacelike coordinate space with 
Euclidean coefficients as the leading term of the metric. 

In Einstein gauge the Weylian scale connection arises from Riemann 

In other 


and alp 




gauge by rescaling with Q = e‘^, (pp = —d^oj (|^ 
words, the additional acceleration due to the scale connection (19) is gener¬ 
ated by the scale invariant representative uj of the scalar field as its potential: 


a^p —Vo; (23) 

If we compare with Newton gravity, we can calculate the hctitious mass den¬ 
sity which one had to assume on the right hand side of the Poisson equation, 
in addition to the real masses, in order to generate the same amount of 
additional acceleration. Obviously here it is 

Pph = {^T^G)-^V'^u. (24) 

In the terminology of the MOND literature the acceleration due to the 
Weylian scale connection corresponds to a phantom energy density We 

see that already on the general level the dynamics of WST differs from Ein¬ 
stein gravity. Only for trivial scalar held, oj = const, the usual Newton limit 
is recovered, otherwise it is modihed. We shall explore how this modihcation 
relates to the usual MOND approaches. 


1.3 WST gravity with cubic kinematic Lagrangian (WST-3L) 

The most common form of the kinetic term for the scalar held is that of 
a Klein-Gordon held, L ^2 = quadratic in the norm of the 

(scale covariant) gradientj^ For our form of the gravitational Lagrangian 
([^ it is conformally coupled for a = —6.^^. Inspired by the approach of the 
relativistic “a-quadratic Lagrangian” (rAQUAL), the hrst relativistic attempt 
of a MOND theory of gravity mm, we hnd that a Weyl geometric scalar 
tensor theory of gravity leads to a MOND-like phenomenology if we add a 
cubic term to the kinetic Lagrangian of the scalar held 


L(j) — L(j)2 + 1^03 • (25) 

The crucial difference to the early approach of rAQUAL is the scale co¬ 
variant reformulation in the framework of Weyl geometry. It results in a 
different behaviour of the scalar held energy density. Bekenstein/Milgrom’s 
model relied crucially on implementing a transition function f{y) between 
the Newton and the deep MOND regime into the kinetic term. The con¬ 
straint of scale invariance of the Lagrangian reduces the underdetermination 
of the Lagrangian and suggests a slightly different form of the kinetic term. 
It is still quite near to the one of the relativistic AQUAL theory. 

^See, e.g., P 48]. 

^In WST Dv denotes the a scale covariant derivative of 0, = dv(j> — 4>'pv 
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In the review paper Bekenstein gives the rAQUAL Lagrangian in the 
form (his eq. (6)) 


= -{87rGN)-^L-^f{L^ S'V'), (26) 


where L is “a constant with dimension of length introduced for dimensional 
consistency” (later it is identified as the MOND acceleration Oq via cP L~^ = 
Oo). Asymptotically f{y) ~ § 2/2 for y <C 1 {MOND regime)] similarly f{y) ~ 
y for y 3> 1 {Newton regime), p is the logarithm of a rescaling function 
between the Jordan frame metric (called the “physical” metric) and the 
Einstein frame (“primitive”) metric y^j/ = Its role is very close 

to our uj in ( |10[ ). A corresponding scale covariant form of the Lagrangian 
could be 

Lct>3 = (^<(')^(?“ V)^ /((C“V)“^ dyUJ d’^uj), (27) 

where in Einstein gauge plays the role of Oq (up to a factor). The sign 

has deliberately been changed; the reasons are given below eq. 


We are here interested in additive modifications (19) of Einstein gravity, 


mainly in a domain in which the effects of the scale connection ipi, = 
compared with those of the Riemannian component of the metric, cannot 
be neglected. This will be called a regime with MOND approximation. For 


non-timelike Vw the Lagrangian (27) becomes 


L4>3 = = (?<(-)'(Jr' <^)"' |Va;|3 

where we have used the abbreviation 

|Va;| := \d^ujd’^u)p 


(28) 


(29) 


(I ... I absolute valu e) t| uj is the scale invariant representative of the scalar 
field introduced in ( |l0| , and V its gradient 

With (po the (constant) value of (p in Einstein gauge we introduce the 
constants 

(30) 


Oio •— ^ (po ) 


On — 


Then the cubic term of the kinetic Lagrangian in Einstein gauge reads 


=-{SttG ao) |Vw|'^. 


(31) 


The dotted equality sign “=”, indicates that the respective equation is not 
scale invariant but presupposes a special gauge made clear by the context 
(similar for f»). Here, as in most cases in this paper, it indicates the Ein¬ 
stein gauge, do has the dimension of inverse length/time and will play a role 

'^For timelike Vo; see |23j . Here we exclusively deal with the spacelike (or null) case. 
^w{4>) = —1 and u)(||Va;|| = —1 imply the scale weight w{L^ 3 ) = —4, as it must be for 
scale invariance of £.^ 3 . 
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analogous to the MOND acceleration Qq ~ [c] H (c the velocity of light). Co¬ 
efficients of type [c] and [h] will often be suppressed in the following general 
considerations. They will be plugged in only in the final step. The pro¬ 
portional factor in Oq = ^ is chosen such that the WST model with cubic 
kinematic Lagrangian (WST-3L) acquires a MOND-like phenomenology in 
a weak gravitational field in which the scalar field and the scale connection 
cannot be neglected. 

is a typical “large number” in the sense of bridging the gap between 
the smallest and the largest physically meaningful energy scales 

^2 ^ ^ ^ io63 _ 

CLo [fi'\ Lpl 

To the kinetic Lagrangian of the scalar field a potential term is added. 
It must be of order 4 to provide for scale invariance of the density (eq. Q: 

Lva = (32) 

The absolute value of the corresponding energy density can be read off from 

in Einstein gauge it is constant. With A = it becomes 

{^4>)~^\Lv4\ = = ^36 0^ = (33) 

comparable to the cosmological constant term A = 3 AIa of the standard 
approach. 

Variation of the Lagrangian leads to the dynamical equations of WST, 
the Einstein equation and the scalar field equation. The scale invariant 
Einstein equation i^ 

Ric-^g = m-^T^^^ + e, (34) 

where g denotes the whole collection of metrical coefficients and the 

energy tensor of matter The scalar field contributes to the total energy 
momentum with two terms, 0 = 0^^^ -|- , the first of which is propor¬ 

tional to the metric (thus formally similar to a vacuum energy tensor) 

0^^^ = <p-^[-DxD^<P^ + r\Lv4 + L^))g (35) 

e™ = r'‘ (36) 

®This means that not only the equation but all its constitutive (additive) terms are 
scale invariant. In particular the Ricci tensor Ric of Weyl geometry is scale invariant 
because the affine connection of the Weyl metric is. 

^ See, e.g. [8], [6l pp. 96ff.]. 
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Varying with regard to (p gives the scalar held equation. Subtracting the 
trace of the Einstein equation for a conformally coupled L ^2 term (a = —6^^) 
strongly simplihes it and introduces the trace of the matter tensor into the 
scalar held equation. In Einstein gauge, with g the Riemannian component 
of the metric, it can be written in terms of the covariant derivative gV with 
regard to g (Levi-Civita connection in Einstein gauge) a^ 

= -47rGaotrr(™) . (37) 

If we introduce the corresponding Riemannian covariant operator 

Pmuj= gVu{\^oj\d''uj) = {d^\Vuj\ d'^u; + \Vu;\ gV^ , (38) 

the scalar held equation for a huid with matter density pm and pressure Pm 
simplihes to the covariant Milgrom equation 


Pm w = dvrG {pm - 3pP ■ 


(39) 


In this derivation, with tr entering by subtracting the trace of the Ein¬ 


stein equation, a sign choice like in (26) leads to the wrong sign on the r.h.s 


of the Milgrom equation. By obvious reasons (38) will be called the covariant 


Milgrom operator. In the static weak held static case cj does not depend on 
the time coordinate. Moreover with g^^^, ~ the expression Vu\d'^p 
turns into the nonlinear Laplace operator V • (| Va;| Vw) of the MOND theory 
with Euclidean scalar product • and norm | ... |. 


In the general case we have to complement (39) with the Einstein equa¬ 
tion in Einstein gauge 


Ric-^g = ST:GT^^^ + e, 


(40) 


In vacuum, the trivial scalar held u = const is a basic solution of (39). 


Then WST reduces to Einstein gravity. In particular, the Schwarzschild and 
the Schwarzschild-de Sitter solutions of Einstein gravity are special (degen¬ 
erate) solutions of WST-3L equations for 4 = 0 or 4 ~ 6, respectively. In 


fact, they solve (40), (39) for cp = const in Riemann gauge, i.e. in the case 


of Einstein gauge equal to Riemann gauge (gp) = (ff, 0). The Riemannian 
component of the metric {g = g =: g) is given by 


ds^ = —(1 - Kr‘^)dt^ -|- (1 - nr"^) ^dr"^ +r‘^{dx\ + siP X 2 dxp^ . 

r r 

(41) 

Then Ric — ^g = —Sng and 0 = 0^^^ = —jpaJ^g. Therefore the Einstein 
equation is satished for 3k = jp dj^, i.e. k ~ 2H^ for ^ ~ 6 and (3 ~ 100. 
We see that in the case of a negligible Weylian scale connection the classical 

^ |23i pp. 15f., sec. 7.2, postprint version arXive v4] 
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(non-homogeneous) point-symmetric solutions of Einstein gravity are valid 
also for the dynamics of WST. This implies that in the case of a negligible 
scale connection Newton dynamics is an effective approximation for point 
symmetric solutions of WST (in Einstein gauge). 

In order to make such a type of Einstein limit compatible with our La- 
grangian, a suppression of the L^a-term for sufficiently large accelerations 
Or of (18) is necessary. Following the example of ordinary MOND theories, 
one might be tempted to plug a factor /(]^^) with a function / such that 

f{y) ~ 1 for y > 0.01 and f{y) ~ 0 for y <C 1 into the r.h.s expression of (28) 
for L (^3 . But such a choice would have the blemish of a coordinate dependent 
argument of the function. A better alternative is provided by the hypothesis 
that the scalar field inhomogeneities are suppressed if any of the sectional 
curvatures k (with respect to the Riemannian component of the metric in 
Einstein gauge) surpasses a certain threshold, e.g., k > (10®ao[c“^])^. In the 
next section we investigate the case of a non-negligible scale connection. A 
more detailed discussion of the transition between the two domains has to 
be left open for another occasion. 


1.4 A WST approach with MOND-like phenomenology 


If the conditions for the weak field approximation (20) are given, it is possible 


to identify a MOND regime as a region in which the Newton acceleration ajv 


is smaller than Oq (here Otv can be identified with or in (22)). Then the 


scalar field equation (39) reduces, in reliable approximation, to 


V • (iVwiVcj) = -IvrGdo , 


(42) 


with the Euclidean V-operator. We call this a MOND approximation. For 
pressure-less matter with energy density pm we get 


V • (|Va;|Va;) = A-nGoo Pm ■ 


(43) 


That is similar to the AQUAE approach]^ Note that only the matter energy 
momentum tensor, without the scalar field energy density, appears on the 


r.h.s. of (42) 


Straight forward verification shows that, independent of symmetry con¬ 
ditions, a solution of (43) is given by uj with a gradient Vcj = —a,p such 
that 


ON 


— \/ho I 


onI 


On 


V |aiv| fiNi 

where ajy denotes the Newton acceleration of the given mass density. 


(44) 


= ArrGpr, 


a]\f = —V4*7v 


(45) 


a, a. 
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(calculations in the approximating Euclidean space with norm | ... |). The 


solution of the non-linear Poisson equation (43) is much simpler than one 


might expect at a hrst glance: In a first step the linear Poisson equation of 
the Newton theory is to be solved, then an algebraic transformation of type 


(44) leads to the acceleration due to the solution of the non-linear partial 
differential equation (43). In fact, has the form of the deep MOND 


acceleration of the ordinary MOND theory (but with a different constant 

do) 


10 


This is only the most immediate modihcation of Newton gravity. In (40) 


there is also the additional term of the energy density due to the scalar field, 
Psf = (87rG)“^0oo- It modihes the r.h.s. of the Newton limit of Einstein 
gravity]^ Neglecting contributions at the order of magnitude of cosmological 
terms (~ H^) and of |Va;p), the energy density of the scalar held in Einstein 
gauge simplihes to 

Psf ^ {A7:G)-^V^uj , (46) 

where Latin indices j,k ... refer to space coordinates only (see appendix 5.1). 

That is equal to the value of the phantom energy density corresponding to 
the acceleration of the scale connection (24). The total “anomalous” additive 


acceleration (in comparison to Newton gravity) is therefore 

^add T ^sf • 

In the central symmetric case 


(47) 


(48) 


In the case of oat <C aadd this leads to MOND-like phenomenology in the 
deep MOND regime if 


Oo = ^ ~ — 77 [c] 2.4 • 10 ® cm s ^ . 


Because of (44) the total acceleration a is then 


a = On + ttadd = Oiv 1 + 


fiNl 




(49) 


(50) 


laivl 


This raises the question of the Newtonian limit. (44) implies |a,^| <C 
in regions where |aAr| Uo (> do)- Therefore can effectively be neglected 
in the case of ‘large’ values of |a 7 v| derived from (|45|). Assuming the hypothe¬ 


sis at the end of section [L3] (or some equivalent), the Newton approximation 


the terminology of the MOND community, the MOND approximation of WST-3L 
behaves like a special case of a QMOND theory [9l pp. 46ff.]. 


^^In contrast psf does not enter the r.h.s. of the scalar field equation (391, and therefore 
does not enter the r.h.s of (45l. 
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is reliable in WST gravity, irrespective of the question of how to characterize 
the transition between the MOND and the Newton approximation. Here 
we shall consider the MOND approximation in an “upper transition” regime 
only, where roughly |aAr| < 10^ Oo 

For centrally symmetric mass distributions p(r) with mass M{r) inte¬ 
grated up to r (where ^ = 1^1 denotes the Euclidean distance from the sym¬ 
metry center, y = {yi,y 2 , Us) the coordinates of the approximating Euclidean 
space) this implies 


12 


= -Vujii - y/GM{r) , 


- 


^/GMl 


r) Uo 


(51) 


Then the phantom energy density (24) become^ 


Pph = (47rG) 


^jGM{r) do 


and also 


Psf ^ (47rG) 


^GM{r) 0.0 


(52) 


1.5 Comparison with usual MOND theories 


We can now compare our approach with other models of the MOND family. 
Simply adding a deep MOND term to the Newton acceleration of a point 
mass is unusual. M. Milgrom rather considered a multiplicative relation 
between the MOND acceleration a and the Newton acceleration oat by a 
kind of ‘dielectric analogy’. 


aN = p {—) a , 
do 


with fi{x) 


or the other way rounc 
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a = z^( —)aAr, with i'{y) 


Oo 


1 

X 

1 

_ 1 
y 2 


for X - 
for X 


for y 
for y - 


oo 
■0 , 


• oo 

0 . 


(53) 


(54) 


From this point of view our acceleration (50) is specified by a well defined 
transition functions 


pLw{x) = 1 + 


1 — -v/TT^ 


2 x 


and Vyj{y) = l + y 2 . 


(55) 


^One might speak of the upper transition regime for Oo < |aiv| < 100 Oo, of the MOND 
regime if |ajv| < flo and of the deep MOND regime for, let us say, |ajv| < 10~^ao |23l sec. 
7.3]. We cannot claim knowledge on the “lower” transition regime with 100 < |ajv| but 

not yet Uo ^ |aiv| (however <C may be specified); see end of section 


1.3 


^^For ds^ = dr^ -|- r^{d6^ + s,iv? 6 dtP) the crucial affine connection components are 
:i = r-h 

X means /i(a:) — x = 0{x), i.e. ^ remains bounded for a: —>■ 0. Cf. 


11 

14 


'‘Here /i(x) 
m 51f.] 


16 


















One has to keep in mind, however, that onr transition functions [i, v are 
reliable only in the MOND regime and the upper transitional regime (ronghly 
On < lO^Oo). They cannot be nsed for discnssing the Newtonian limitd^ It 
will be important to see how they behave in the light of empirical data, in 
particnlar galactic rotation cnrves and clnster dynamics. 

In the MOND literatnre the amonnt of a (hypothetical) mass which in 
Newton dynamics wonld prodnce the same effects as the respective MOND 
correction Oadd is called phantom mass Mp^. For any member of the MOND 
family the additional acceleration can be expressed by the modified transition 


fnnction u = u — 1 with u like in (54) 


O^add — 


|aiv| 


ON ■ 


(56) 


The phantom mass density pph attribnted to the the potential ^ph satisfies 
AttG Pph = S/^^ph and V4>p/j = —Oadd- A short calcnlation shows that it 
may be expressed as 


Pph 


= V 


Qjvl 

CLn 


Pm - (dvrGao) 


|QAf| 

(Xo 


(V|aAr|) • OAT. (57) 


It consists of a contribntion proportional to pm with factor u, which dom¬ 
inates in regions of ordinary matter, and a term derived from the gradient 
of |a 7 v| dominating in the “vacnnm” (where however scalar field energy is 
present). For the Weyl geometric model with Pw{y) = y~^, P'wiy) — ~\y~^ 
Pph tnrns into: 

Pt 
Psf 


|aiv| 


Pm + (SvrG) 


-b 


Vl«Af| 


V(|a^|)- 


|0Ar| 


Pph — 2 P^ 


(58) 

(59) 


The first expression of (59) is compatible with (46). 


In onr case it wonld be ntterly wrong to consider the whole of pt as 
“phantom energy”. Half of it are dne to the scalar field energy density, 
the scalar field halo pgf, and expresses a true energy density. This energy 
density appears on the right hand side of the Einstein eqnation (40) and the 


Newtonian Poisson eqnation as its weak field, static limit. It is decisive for 
lensing effects of the additional acceleration. The other half, pph, is phantom, 
i.e. a fictitions mass density prodncing the same acceleration as the Weylian 


scale connection (24). Only for the sake of comparison with other MOND 


models we may speak of pt as some kind of gross phantom energy,in contrast 
to the “net” phantom energy pph i. 


^®See fn. 12 and the text above it. 
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We have to distinguish between the influence of the additional structure, 
scalar field and scale connection, on light rays and on (low velocity) trajec¬ 
tories of mass particles. Bending of light rays is influenced by the scalar field 
halo only, the acceleration of massive particles with velocities far below c by 
the scalar field halo and the scale connectionF^ 

Also in another respect our theory differs from the usual MOND ap¬ 
proaches. In MOND external acceleration fields of a system under consid¬ 
eration are difficult to handle. In WST, like in GR, a freely falling (small) 
system does not feel the external acceleration field if it is sufficiently small, 
relative to the inhomogeneities of the external gravitational field, for ne¬ 
glecting tidal forces. In this sense, the external acceleration problem does 
not arise in the WST MOND approximation (42). 

Another important consequence follows: the scalar field energy formed 
around a freely falling subsystem of a larger gravitating system, calculated in 
the MOND approximation of the freely falling subsystem, contributes to the 
r.h.s. of the Einstein equation of any other subsystem (in relative motion) 
and also to that of a superordinate larger systemj^ This has to be taken 
into account for modelling the dynamics of clusters of galaxies. 


1.6 Short resume 


We have derived the most salient features of the Weyl geometric MOND 
approximation (WST MOND) and are prepared for a comparison with em¬ 
pirical data. Before we do so, it may be worthwhile to collect the results 
which are necessary for applying it to real constellations in a short survey. 

Consider a gravitating system which in the Newton approximation of 
Einstein gravity is described by the baryonic matter density pm, the accel¬ 
eration Um and potential with 

= IvrGpm, am = -V^m- (60) 


The modification due to WST MOND leads to an additional acceleration 
Oadd with the following features: 

(i) The total acceleration a is a = Om + aadd with 


a = a 



= V 


l^m I 

Q/n 


(61) 


where v{y) = 1 + y 2 . For Om 3> ao the Newton approximation 
applies. (61) holds for |am| < lO^Oo only (“upper” transition regime). 


^®That may look like bad news for explaining lensing at clusters and microlensing at 
substructures. But the particular transition function seems to compensate much of this 

effect. _ 

^^In principle that presupposes that the whole energy momentum tensor (35 361 (and 
its system dependent representation) is considered. For slow motions and weak field ap¬ 
proximation a superposition of energy densities like in Newton dynamics seems legitimate. 
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No information can be drawn from it for \am\ larger bnt not yet S> Oo 
(the “lower” transition regime). 

(ii) The “reciprocal” transformation fnnction defined by am = (|a^) ® 

is 

/ N 1,1“ Vl + , , 

Kx) = 1 +- ^ -■ (62) 

(iii) aadd consists of two components Uadd = + Os/ = 2o;^. The first one 

is derived from a potential w satisfying the non-linear Poisson eqnation 

V • (I VculVa;) = ttGoo Pm ) 0^, = —Vw. (63) 


(iv) The second one, agf, can be nnderstood as a Newton acceleration dne 
to the energy density psf of a scalar field (part of the modified gravita¬ 
tional strnctnre). Its potential satisfies a Newtonian Poisson eqnation. 
It satisfies 

- Va*/ = ^-nGpsf (64) 

with energy density 

fpm + (87rG)"^V(|am|) • 7^) (65) 

2 \\am\J V l«ml/ 

Psf is part of the energy-momentnm tensor of the scalar field cj) and in 
this sense “real” rather than phantom. 


(v) 


is formally derivable in Newton dynamics from a fictitions energy 
density 

Pph — Psf ■ (66) 


Pph is the net phantom energy of WST MOND. For comparison with 
other models of the MOND family one may like to consider pph + Psf = 
pt as a kind of “gross phantom energy” (althongh the larger part of it 
is real). It is transparent rather than “dark” (see (74) below). 


(vi) (i) - (v) are reliable approximations also for small (local) gravitating 
systems freely falling in a larger gravitating system, if am ~ 10 ^ Oq in 
the local system. The snbsystem can be considered as “small” with 
respect to the snper-system, if tidal forces of the snper-system can 
be neglected. In hierarchical systems like galaxy clnsters the energy 
density contribntions psf of the snbsystems and the snper-system (cal- 
cnlated in different reference coordinate systems) add np to the total 
energy density of the scalar field, if the velocities of the snbsystems rel¬ 
ative to the barycenter of the snper-system are small. This is a crncial 
difference between WST MOND and ordinary MOND theories. If one 
likes, Psf can be considered as the “dark matter” component of WST 
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MOND; although as the energy density of the scalar held it is not 
constituted by the usual (hypothetical) quantum particles (WIMPs, 
axions etc.). To demarcate this difference it might better be called 
transparent matter/energy of WST. 

(vii) Gravitational leasing is due to the scalar held energy density only {psf), 
while the dynamics of WST corresponds to the total phantom density 
{pt = Psf + Pph)- It remains to be seen whether such a difference is in 
agreement with observations. 


2 Halo model for clusters of galaxies 

2.1 Cluster models for baryonic mass (hot gas and stars) 

In the astronomical literature, the density prohle of hot gas and (smeared) 
star/galaxy matter in a galaxy cluster is often described by a centrally sym¬ 
metric prohle of the following form: 


p{r) = po 




(67) 


/3 is the ratio of the specihc energies of the galaxies and the gas, po the central 
density and Tc is the core radius da 1201 [igR ( |67[ ) is called a fi-model for 
the mass distribution. For our test we assume density models for the gas 
mass Pgas{f) and for the galaxy mass Pstar{f) with the same form parameters 
/3 and rc. We thus work with an idealized model using proportional density 
prohles for the hot gas and for the galaxies with parameters (3, Vc determined 
from observations of the hot gas. The central densities po can, in principle, be 
determined from mass data for gas Mgasiri), respectively stars Mstar{fi)^ at 
a given distance ri. The empirical determination of Mgas{ri) and Mstari^i) 
from directly observable quantities is a subtle question; it will be discussed 
in section l3T] 

Large scale gravitational effects on the cluster level are often modelled 
in the Newton approximation of Einstein gravity with baryonic matter and 
an assumed dark matter halo which is inferred from its gravitational effects 
(in Einstein/Newton gravity). Another, minority, approach in the literature 
works with an evaluation of the data in a MOND limit of the most well known 
relativistic MOND theory TeVeS. R. Sanders is one of its protagonists; he 
concludes that in this approach a much smaller amount of unseen matter has 
to be assumed in addition to the baryonic mass. Its value is consistent with 
the hypothesis of a halo of sterile neutrinos, concentrated about the cluster 
center [JU]. Here we want to explore the feasibility of the WST approach, in 


is the distance from the cluster center at which the projected galaxy density is half 
the central density po- 
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particular regarding the question of how much unseen matter has to be added 
to the gravitational effects of the model in order to reproduce (“predict”) the 
observed accelerations, respectively their measurable effects. 


2.2 Two contributions to the scalar field energy in clusters of galaxies 


The mass distribution of the hot gas pgas in a galaxy cluster is described 
by a /3-profile (67) in a locally static coordinate system with origin at the 
barycenter of the cluster. The averaged star mass will be represented by a 
continuous distribution Pstar of a /3-profile with the same parameters (3, Vc, 
but with a different value of po- Both together form a continuity model of 
the baryonic mass distribution pbar = Pgas + Pstar- Estimates show that 
Newtonian gravitational accelerations induced by pbar (fnr away from mass 
concentrations stars, galaxies, and galactic centers) are below 10^ • Qq- They 
are small enough for allowing to in a weak field static approximation with 


the scalar field equation in the MOND approximation (43). The resulting 
contribution to the scalar field energy will be called Psfi- The additional 


acceleration of the Weyl geometrical scale connection (66) can be expressed 


in terms of a phantom mass density which will be called pph i. 

In this first approximation the star mass is approximated on a par with 
the hot gas, i.e., it is described by its continuously smeared out mean den¬ 
sity. But stars are agglomerated in galaxies which form freely falling subsys¬ 
tems of the cluster with considerable inter-spaces in the super-system (the 
cluster). For each subsystem a locally static coordinate system with origin 
at the respective galactic center can be chosen. In this system the local 
inhomogeneities of star mass distribution in the cluster, and the resulting 
inhomogeneities of the gravitational field, in the vicinity of of the galaxy can 
be calculated. On the galaxy level the MOND theory has proven effective for 
modelling gravitational effects deviating from Einstein and Newton gravity 
without assuming real dark matter [9] . Although we expect that the MOND 
approximation of WST gravity shows similar features, this is not the point 
in the present investigation. 

Here we are interested in the neighbouring regions of galaxies as sub¬ 
systems of their respective cluster. These subsystems form scalar field halos 
of their own which contain real energy (different from the classical MOND 
theory which leads to phantom halos only). In the framework of the present 
approach, the scalar field halos in the neighbourhood of each galaxy con¬ 
tribute to the energy density which adds up globally, i.e. on the cluster 
level, to a component of scalar field energy psf 2 which has been suppressed 
in the first continuity approximation of the total baryonic mass. In a second 
step we therefore determine this component approximately and add it to the 
total the scalar field halo. 
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2.3 Scalar field and phantom halos Psfi, Pphi in the cluster-barycentric 
MOND approximation 

The baryonic mass up to radius r, 


Mbar{r) = dvr / Pbar{u) V?du , 

Jo 


( 68 ) 


determines the Newton acceleration Ubar = due to the total baryonic 

mass. The densities of the scalar field halo psfi and the phantom halo of 
WST MOND follow from (65), (66). They are 


Psfi = 


1 


2 V I dbar I 
Pph 1 = Psfi ■ 


Pm + iSnC) V(|a6ar|) 


^bar 

I ^har I 


(69) 

(70) 


The respective masses of the halos Mgfi, Mphi arise from integration. 


2.4 


Scalar field halos of galaxies in their respective galacto-centric 
MOND approximations 


As already indicated in section 2.2 (69) and (70) do not make allowance 


for the fact that the star matter forms a discrete structure of an ensemble 
of galaxies each of which is falling freely in the inertio-gravitational field of 
the super-system (hot gas and other galaxies). Every galaxy possesses a 
local MOND approximation with regard to its own barycentric static ref¬ 
erence system. The acceleration Obar of the super-system (with respect to 
the barycenter rest system of the hot gas) is transformed away in each of 
the local MOND approximations. The latter leads to a galactic scalar field 
halo which persists under changes of reference systems with small, i.e. non- 
relativistic, relative velocities. It contributes to the total energy of the scalar 
field, calculated in the cluster barycentric system. (Of course this is not the 
case for the phantom halo of the single galaxies.) In principle, we have to add 
up all these effects to a scalar field energy density psf 2 in order to fill in this 
lacunae. But an exact calculation would have to solve a highly non-trivial 
A^-body problem for the motion of the galaxies. 

The experience with the calculation of the combined MOND halos of stars 
inside galaxies shows that a resolution of the star matter inside galaxies 
down to individual galaxies is not necessary to achieve good results. In 
the outer region of galaxies a continuity model for the distribution of star 
matter in the galactic disk gives reliable approximations for the MOND 
acceleration. Similarly we want to check whether also here a continuity 
model for the system of galaxies alone, abstracting from the gas mass, leads 
to an acceptable approximation for Psf 2 - For this calculation, the gas mass 
has to be omitted because the galaxies are falling freely in the outer field of 


22 











the cluster; the gravitational potential of the latter does not enter the local 
MOND approximation of the galaxies. 

Using (65) again we get for the second (inhomogeneity) component of 
the scalar held energy 


Psf2 = 


1 


O^star 


Pstar (SttG) V(|flstar|) 


(^star 

I Oistar I 


(71) 


with |astar(^)| = and Mstar{r) the integral analogous to (68) for 

the star density pstar- 

We hnally arrive at a halo model for galaxy clusters constituted by the 
components psf i, Psf 2 , Pph i ■ All of them are determined by the two compo¬ 
nent baryonic prohle of the cluster. 


2.5 A three-component halo model for clusters of galaxies 


In addition to the Newtonian gravitational effects of the baryonic mass den¬ 
sity 

Pbar — Pgas Pstar j C^^) 


modeled by a /3-model of type (67), the WST MOND approach predicts ac¬ 


celerations generated by the scalar field halos (@, (0. Because of their 
small values their combined effect can be approximated by a linear superpo¬ 
sition in the barycentric reference system of the cluster]^ 


Psf ^ Psf 1 + Psf 2 ■ 


(73) 


Moreover, there arises an acceleration due to the scale connection in 


the barycentric rest system of the cluster (45). Its gravitational effects are 


representable by the hctitious (net) phantom halo Pphi of (70) 


Pph i — Psf i- 

On the other hand, the phantom energies of the individual freely falling 
galaxies do not survive the transformation to the cluster rest system and do 
not play a role on the cluster level. 

In the usual MOND theories there is no scalar held energy; all additional 
effects with regard to Newton dynamics may be ascribed to a (hctitious) 
phantom energy density. Phantom energy densities of single galaxies do not 
survive the transformation to the cluster barycentric system. In MOND there 
is therefore no analogy to Psf 2 ] the latter is the crucial distinctive feature 
between the two approaches. For a comparison of WST-3L and usual MOND 

^®Because of slow (i.e., non-relativistic) relative velocities of the galaxies the energy 
densities of their respective scalar field halos can be taken over to the cluster barycentric 
reference system. 
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approaches with regard to galaxy clusters it is not sufficient to evaluate the 
difference between the transformation functions ffix) (55). 

For an even wider comparison with other approaches it may be useful to 
add up the scalar held and phantom halos to a kind of “dark matter” halo or, 
more precisely, to a substitute for the latter. But one must not forget that 
in the present model there is no dark matter in the ordinary sense. Here we 
only hnd a transparent halo made up of the (real) energy density of the scalar 
held and the (hctitious) phantom energy density ascribed to the acceleration 
effects of the scale connection in Einstein gauge (with respect to the cluster 
barycentric rest system): 

Pt = Psf + Pphi (74) 


From the gravitational leasing point of view, it would be even more appro¬ 
priate to consider psf alone as the WST equivalent of a dark matter halo, 
not forgetting that even the real halo Psf is not due to fermionic particles, 
but to the scalar held, and thus to the extended gravitational structure of 
WST. From a quantum point of view, the scalar held has to be quantized if 
one wants to search for a (bosonic) particle content of pt or psf. 

The total dynamical mass of the model (up to some distance r from the 
center of the cluster) is 


Mtot — Mhar + Mt, 


Mt = Msf + Mph 1, Msf = Msf 1 -I- Msf 2 


(75) 


with Mfyar = Mgas + Mgtar- The lensing mass is a bit smaller. 


Miens = Mbar + Mgf . 


(76) 


Mathematically, the integral of the scalar held energy density to arbitrary 
distances diverges. Like in the dark matter approach a virial radius of a 
cluster may be dehned, which roughly delimits the gravitational binding 
zone of the cluster. Comparing pt (respectively {psf)) with the critical energy 
density pcrit of the universe, one may, for example choose r 2 oo with 


Ptij'20o) ~ 200 Pcrit 


(77) 


as a representative of the virial radius. 

Not far beyond the gravitational binding zone of the cluster, the energy 
density will have fallen to such a small amount that its centrally symmet¬ 
ric component is inconceivably stronger than the density huctuations in the 
inter-cluster space. To continue the integration into this region, and beyond, 
has no physical meaning. In the long range the energy density of the scalar 
held approaches the cosmic mean energy value. A physical limit of integra¬ 
tion has to be chosen close to the virial radius beyond which the gravitational 
binding structure of the cluster is fading out. 
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3 A first comparison with empirical data 

3.1 Empirical determination of mass data for galaxy clusters 

For a first empirical exploration we confront the WST cluster model with 
recent mass data for 19 clusters obtained on the background of Einstein/New¬ 
ton gravity and ACDM by a group of astronomers about Y.-Y. Zhang, T.F. 
Lagana and T. Reiprich EH, m- We use the form parameters /3,rc of the 
/3-models of these clusters, published in an earlier study by one author of 
the group [18]. In the present study we have to take the mass data and the 
form parameters essentially at face value. Methodological questions arising 
from this procedure are discussed in the next subsection. There seem to be 
sufficient reasons for expecting that the different background theories do not 
principally invalidate the results thus obtained. Of course, an authoritative 
empirical study would presuppose an evaluation of observational data on the 
background of WST itself; it can be done only by astronomers, if they get 
interested in the present approach. 

The studies of Zhang/Lagana/Reiprich e.a. have the great advantage to 
build upon three independent observational data sets for determining the gas 
mass Mgas,500) the star mass M^^soo and the total mass M 500 (assuming a 
dark matter explanation for the observed gravitational effects) at the refer¬ 
ence distance rsoo- The latter is determined for each cluster at the distance 
rsoo from the cluster center at which the total gravitational acceleration 
indicates a total mass density 500 times the critical density. 

(i) Mgas,500 has been extracted from X-ray data on the hot intraclus¬ 
ter medium (ICM) collected by XMM Newton and ROSAT. Surface 
brightness data have been used to infer an ICM radial electron num¬ 
ber density profile, and spectral analysis data gave information on the 
radial temperature distribution. From that a gas density distribution 
has been reconstructed and the gas masses at rsoo by integrationj^ 

(ii) M*^ 5 oo has been determined from optical imaging data due to SDSS 7 in 
two steps. First the total luminosity of the cluster has been determined 
by means of a “galaxy luminosity function” (GLF); then the mass is 
estimated using mass-to-light ratios depending on the cluster mass. In 
the last step models of the star development in the respective galaxy, 
elliptical or spiral, enter. They depend on assumptions on an “initial 
mass function” (IMF). Two possibilities for the IMF (Salpeter versus 
Kroupa) are considered and compared in [27j . |28) . According to the 
authors the difference of the stellar mass estimate can result in a factor 

2 m p-4]. 

(iii) The total cluster mass M 500 has been determined on the basis of the 
velocity dispersion of galaxies, using spectroscopic data from |26[ tab. 

outline of the procedure and literature for more details is given in |27l 3]. 
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1]. The mass estimator used is eq. (2) of [S] 

A = 1.50 ± 0.02, (7y the 3-dimensional velocity dispersion inside a 
sphere of virial radius (by convention = r 2 oo)- Reasons for this 
choice are given in |5l sec. 3]. M 500 was then determined from M 200 
by a NFW-model. 

3.2 Theory dependence of mass data for galaxy clusters 

Mass densities of the hot gas and of star matter in galaxy clusters are in¬ 
directly inferred from observable quantities; they are thus theory dependent. 
Even inside the same background theory they may depend on choices of 
models and methods of evaluation. 

That makes it a difficult task to compare our model with empirical data. 
A fine-grained judgement presupposes an evaluation of observational raw 
data on the background of WST gravity or, at least, a detailed estimation 
of systematic errors resulting from a comparison of different background 
theories (Einstein gravity with ACDM and Newton approximation, or al¬ 
ternatively TeVeS-MOND, in comparison with WST and its MOND approx¬ 
imation). This task has to be left to astronomers, if they become sufficiently 
interested in the present approach. But, taking this caveat in mind, it still 
seems possible to confront available data from, e.g., the Einstein gravity- 
ACDM-Newton approximation framework with our model, in order to get 
a first impression of its potential usefulness. A comparison with mass data 
derived in a TeVeS-MOND background would give welcome supplementary 
information. This is not attempted here. 

Let us discuss the possibility and the problems of a confrontation of these 
data with the MOND-approximation of WST: 

(1) The mass of the hot gas (intracluster medium) Mgas 500 (at ^ 500 ) has 
been determined in the mentioned study from X-ray data obtained by 
XMM-Newton and ROSAT. The temperature of the gas is estimated by 
a fit to the measured spectrum. The gas density p{r) is reconstructed, 
using model assumptions, from intensity observables and then inte¬ 
grated up to r5oo0 Up to usual model dependence, the transfer of the 
mass data from the standard gravity background to WST seems to be 
relatively uncritical. 

(2) Several methods for determining the stellar mass M^soo are mentioned 
in |28]. In this study the star mass is gained from optical imaging data 
due to SDSS 7 in two steps indicated in (ii) above. According to the 

this evaluation the hydrostatic assumption was corrected by taking the velocity 
dispersion into account |28l sec. 2.2f.]. 
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authors the difference of the stellar mass estimate due to different initial 
mass functions can result in a factor 2 m p.4 ]. Another approach 
would be to estimate stellar masses of the individual galaxies and “to 
construct the stellar mass functions in order to sum the stellar masses” 
(ibid, 1). Moreover, an additional component of star matter can be 
associated to the intracluster light. All in all, the estimate of the 
star mass concentrated in galaxies seems to depend more on models of 
galactic star evolution than on the background gravity theory. In spite 
of that the precision cannot be expected to be better than by a factor 
2 (respectively 0.5). 


(3) In Einstein gravity / AC DM the cluster mass can, in principle, be esti¬ 
mated from the velocity dispersion a of galaxies at distance R (from the 
center) by an estimator derived from the virial theorem M ~ G~^a^R. 


The additional acceleration of WST-3L aadd = a^p + agf (47) is dynami¬ 
cally indistinguishable from the effects of “true” Newtonian masses. So 
far it seems as if the estimation of total mass can be transferred to the 
MOND approximation of WST without problems. But if the radius 
R does not include the “whole” cluster mass (however defined), as is 
here the case (item (iii) last subsection) a surface pressure term must 
be taken into account. That complicates the case. 

In standard gravity the necessary correction is implemented by a cubic 


mass estimator M ~ given above as (78). Moroever, the authors of 


our reference study m reconstruct M 500 and rsoo from these values 
using the NEW profile. Because of the different profile for the scalar 
field halo of WST this is a critical step for our exercise On the 


other hand, if the resulting systematic errors are smaller than the error 


intervals of Mtot (75), implied by the observational errors of the other 


quantities, they do not disturb a rough empirical check of the model. 


(4) Finally the dependency of the data evaluation on the background cos¬ 
mology has to be taken into account. The data of the 19 clusters used 
in the following have redshift z < 0.1. The geometrical and dynami¬ 
cal corrections implied by the AC DM cosmology are correspondingly 
small. An evaluation in, e.g., a Lemaitre-de Sitter model (or even a 
non-expanding Weyl geometric model with redshift) |231 section 4.2] 
would affect the data only by a minor expansion of the error intervals. 


The points (1), (2) and (4), in particular the estimate of stellar mass concen¬ 
trated in galaxies and gas mass, are fairly insensitive against a change of the 
background theory from Einstein gravity to WST. The theory dependence 
of rsoo is uncritical in our context. Any other reference radius could have 

ex-post comparison of the NFW-halo and the WST-halo for the Coma cluster is 
given in fig. 
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been taken, as long as it is specified in astronomical distance units. The 
estimate for the total masses at rsoo and r 2 oo is the critical point for our 
purpose (item (3)). However, if the difference of the halo profiles between 
the scalar field energy density of WST and NFW dark matter does not push 
the estimates for the total masses M 500 , ^^200 outside the error intervals of 
our halo model (due to observational input data), we may still be able to 
draw first inferences from the following evaluation. 

3.3 Empirical data for 19 clusters 

The studies m, m contain new data on the baryon content and the 
total gravitational mass for 19 clusters of galaxies (as it appears in an 
Einstein-ACHM framework with Newton approximation) The mass data 
M 500 , Afgassoo, Af* 500 are given in columns (5), (7), (8) of table 1 in |2^. It 
is reproduced in our fig. The values for r 5 oo are published in |25l tab. 1, 
col. (5)] (here table 1). A comparison of the total cluster masses derived 
from the velocity dispersion with a mass estimate derived from the gas mass 
shows that the two clusters A2029 and A2065 are outliers, with total clus¬ 
ter masses considerably higher than the corresponding gas masses would let 
expect. The authors therefore separate the two outliers from the rest of the 
data, with the remaining 17 clusters as a reliable data set [271 P- 3]. We 
shall do the same. 


Name 

X-ray center (J2000) 

RA Dec 

Redshift 

Msoo 

10^^ Mo 

M500,M-Argas 
10^‘^ Mo 

Mgas.500 

10^^ Mo 

M*,500 

10*2 

Undisturbed 
/cool core 

A0085 

00:41:50.306 

-09:18:11.11 

0.0556 

6.37 ± 1.00 

5.68 ±0.37 

8.13 + 0.38 

7.36 ± 1.00 

Y/S 

A0400 

02:57:41.349 

+06:01:36.93 

0.0240 

1.83 ±0.39 

1.07 ±0.07 

1.36 ±0.05 

4.39 ± 1.06 

N/N 

lIlZw54 

03:41:18.729 

+ 15:24:13.91 

0.0311 

1.91 ±0.58 

1.18 ±0.08 

1.45 ±0.26 

4.57 + 0.56 

Y/W 

A1367 

11:44:44.501 

+ 19:43:55.82 

0.0216 

1.76 ±0.27 

2.11 ±0.14 

2.07 ± 0.07 

4.35 ± 0.74 

N/N 

MKW4 

12:04:27.660 

+01:53:41.50 

0.0200 

0.50 ±0.14 

0.58 ± 0.04 

0.47 ± 0.02 

1.16±0.22 

Y/S 

ZwC11215 

12:17:40.637 

+03:39:29.66 

0.0750 

4.93 ± 0.98 

4.34 ±0.28 

6.10 ±0.29 

7.05 + 0.83 

Y/N 

A1650 

12:58:41.885 

-01:45:32.91 

0.0845 

3.44 ± 0.66 

4.28 ±0.27 

5.09 ± 0.73 

7.47 ± 1.13 

Y/W 

Coma 

12:59:45.341 

+27:57:05.63 

0.0232 

6.55 ±0.79 

6.21 ±0.40 

8.42 ± 0.63 

13.14 ± 1.80 

N/N 

A1795 

13:48:52.790 

+26:35:34.36 

0.0616 

3.41 ±0.63 

4.46 ±0.29 

5.11 ±0.14 

6.21 ±0.98 

Y/S 

MKW8 

14:40:42.150 

+03:28:17.87 

0.0270 

0.62 ±0.12 

1.10 ±0.07 

0.80 ±0.12 

1.61 ±0.23 

N/N 

A2029 

15:10:55.990 

+05:44:33.64 

0.0767 

14.70 ±2.61 

6.82 ±0.44 

13.35 ± 0.53 

9.59+ 1.11 

Y/S 

A2052 

15:16:44.411 

+07:01:12.57 

0.0348 

1.39 ±0.28 

2.03 ±0.13 

1.86 ±0.10 

3.53 ±0.40 

Y/S 

MKW3S 

15:21:50.277 

+07:42:11.77 

0.0450 

1.45 ±0.34 

2.29 ±0.15 

2.13 ±0.09 

3.90 ± 0.43 

Y/S 

A2065 

15:22:29.082 

+27:43:14.39 

0.0721 

11.18± 1.78 

3.35 ±0.22 

7.66 ± 1.44 

7.32 ±0.75 

N/W 

A2142 

15:58:19.776 

+27:14:00.96 

0.0899 

7.36 ± 1.25 

10.26 ±0.66 

13.76 + 0.73 

8.42 ± 0.77 

Y/W 

A2147 

16:02:16.305 

+ 15:58:18.46 

0.0351 

4.44 ± 0.67 

3.63 ±0.23 

5.04 ± 0.53 

6.84 + 0.90 

N/N 

A2199 

16:28:37.126 

+39:32:53.29 

0.0302 

2.69 ± 0.42 

2.64 ±0.17 

2.97 ± 0.30 

4.76 ±0.50 

Y/S 

A2255 

17:12:54.538 

+64:03:51.46 

0.0800 

7.13 ± 1.38 

4.08 ± 0.26 

7.11 ±0.33 

6.74 ± 0.97 

N/N 

A2589 

23:23:56.772 

+ 16:46:33.19 

0.0416 

3.03 ±0.75 

1.88 ±0.12 

2.54 ±0.17 

5.12±0.56 

Y/W 


Notes. The cluster mass, M5oo,M-A/gas ^ is derived from the Msoo - Mgas,5oo relation, and only used for comparison with the cluster mass, M500, derived 
from the “harmonic” velocity dispersion. “S”, “W”, and “N” denote strong cool-core, weak cool-core, and noncool-core clusters. 

Figure 1: Tab.l of |28) . Properties of the 19 galaxy clusters 


^^|28) contains a correction to the first paper. Here, of courese, we use the corrected 
data. 


28 







Parameters (/3, rc) for the models of these galaxy clusters (as well as of 
many more) have been published earlier in |18[ tab. 4.1], This publication 
also contains mass data M 200 , Mgas 200 at ^200 (with error intervals), mass 
values Ma, MgasA (without error interval) at the Abell radius, here defined 
as TA = 2.14 Mpc, but no data for star masses ^ In |18) the methods for 


determining the total mass and the gas mass were not yet as refined as in 
the later study. It is therefore not possible to aggregate the different data 
sets to one coherent ensemble, As no updated values for the parameters 


(/3, Tc) of the mass profiles were available to the present author, the values 
of |18| are here used as estimators for the form parameters of the /3 model. 

We consider M 500 , Mgas 500 , M^sqo, from |27l|28] as our crucial mass 
data (including reference radii). Mass data at r 2 oo (as well as r 2 oo itself) from 
the older study are welcome as additional information] but they will not be 
used as core criteria for the empirical test of our model. Table 1 collects the 
data used, It shows that the cluster ensemble covers an order of magnitude 


variation for the gas mass Mgasboo and one and a half orders of magnitude 
variation in total mass The selection method by intersecting the clus¬ 

ter sets of different raw data sources does not seem to be influenced by any 
particular bias. We thus may consider the collection as a reasonable data 
set for testing our cluster halo model. 


3.4 Adapting the WST halo model to the data 


For the construction of the model we have to work in Einstein gauge (|^. 
Consistency with the deep MOND acceleration demands ( |49[ ). Taken to¬ 
gether these conditions fix the coefficients of the Lagrangian Q and (28), 
independently of the convention chosen for (fo- With a value for A on the 


order of magnitude 10 (e.g., 4=6 like at the end of section 1.3) the con¬ 


tribution of the Ly 4 term lies many orders of magnitude below the energy 
densities the dominant term in ( p^ and is negligible in our context. 

In agreement with section]^ the test of our halo model for each of the 19 
galaxy clusters can now proceed as follows: 


(I) Specification of the /3 models (67) for gas and for star mass; parameters 
(/3, Tc) from [18], Po determined by fitting to Mgas 500 and M* 5 oo at rsoo 
respectively l28l 


Evaluated for the value of Ho assumed in the later publications |27) . [28| . h = 0.7. 

^®The values for M 500 and Mgassoo given here differ from the ones in outside the 
error intervals. 

^®For error intervals see the respective source, [28| . m- 

^^If one wants to investigate the external halo (beyond r 2 oo) one has to choose a cutoff at 
r 2 oo or fading out functions. The main part of of our investigation deals with the internal 
halo; only in the final discussion, section]^ questions of the external halo come into the 
play. This remark applies, mutatis mutandis, to items (3), (4) below 
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Table 1. Data set used for halo model (error intervals omitted) 


Cluster 

/5 

Tc 

^500 

M 500 

Mgas 500 

A7*5oo 

M 200 

^^200 

Ma 

Coma 

0.654 

246 

1.278 

6.55 

8.42 

13.14 

13.84 

2.3 

12.86 

A85 

0.532 

59.3 

1.216 

6.37 

8.13 

7.36 

7.71 

1.9 

8.72 

A400 

0.534 

no. 

0.712 

1.83 

1.36 

4.39 

1.48 

1.09 

2.93 

IIIZw54 

0.887 

206 

0.731 

1.91 

1.45 

4.57 

2.81 

1.35 

4.51 

A1367 

0.695 

274 

0.893 

1.76 

2.07 

4.35 

4.06 

1.53 

5.77 

MKW4 

0.44 

7.86 

0.58 

0.5 

0.47 

1.16 

0.71 

0.86 

1.79 

ZwC1215 

0.819 

308 

1.098 

4.93 

6.1 

7.05 

10.37 

2.09 

10.65 

A1650 

0.704 

201 

1.087 

3.44 

5.09 

7.47 

11.14 

2.15 

11.11 

A1795 

0.596 

55.7 

1.118 

3.41 

5.11 

6.21 

10.99 

2.14 

11.04 

MKW8 

0.511 

76.4 

0.715 

0.62 

0.8 

1.61 

2.38 

1.28 

4.0 

A2029 

0.582 

59.3 

1.275 

14.7 

13.35 

9.59 

13.42 

2.29 

12.59 

A2052 

0.526 

26.4 

0.875 

1.39 

1.86 

3.53 

2.21 

1.25 

3.79 

MKW3S 

0.581 

47 

0.905 

1.45 

2.13 

3.9 

3.46 

1.45 

5.11 

A2065 

1.162 

493 

1.008 

11.18 

7.66 

7.32 

16.69 

2.45 

14.46 

A2142 

0.591 

no 

1.449 

7.36 

13.76 

8.42 

15.03 

2.36 

13.61 

A2147 

0.444 

170. 

1.064 

4.44 

5.04 

6.84 

3.46 

1.45 

5.15 

A2199 

0.655 

99.2 

0.957 

2.69 

2.97 

4.76 

4.80 

1.62 

6.37 

A2255 

0.797 

423 

1.072 

7.13 

7.11 

6.74 

13.32 

2.27 

12.52 

A2589 

0.596 

84.3 

0.848 

3.03 

2.54 

5.12 

3.58 

1.47 

5.24 


Source: first and last block |18| . middle block |28) . rsoo |25j 

Units: Vc in kpc, r 5 oo,r 2 oo in Mpc, M 500 in 10^"^ Mq, Mgasboo in 10^^ Mq, 

iff*500 in 10^^ Mq, 
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(2) Determination of the Newton accelerations of the baryonic mass com¬ 
ponents. 


(3) Calcnlation of the scalar field halo and the phantom halo of the bary¬ 
onic mass by (69), ( 66 ). 

(4) Calcnlation of the scalar field halo of the system of freely falling galaxies 


(5) Aggregation of these to the total halo (74); choice of fading ont fnnc- 
tions beyond r 2 oo (see appendix). 

( 6 ) Integration of the densitities to the corresponding masses: scalar field 
energy of the galaxies Msf 2 , of the gas mass Mgfi, total scalar field 
halo Mgf, net phantom energy of the baryonic mass Mphi, finally the 
total transparent matter Mf and the lensing mass Mi^^ns (75), (76). 

(7) Calcnlation of the error intervals of the model at selected distances 
(^’500, ?’20o)- 


( 8 ) Comparison of the empirical valne for M 500 (respectively M 200 ) with 
the model valne Mt(r 5 oo) (resp. Mt(r 2 oo))- 

The resnlt of this test are given in snbsection |3.6| Before we tnrn to the 
overall evalnation, we shall have a look at one clnster as an exemplary case. 


3.5 The WST halo model with the Coma cluster as test case 

For a first check of onr model we choose the Coma cluster. According to item 
( 1 ) of the last section we nse the empirical inpnt data from table 1 (nnits 
given there): 

i/3,r„r5oo,Mgas500,M,50o) = (0.654, 246, 1.278, 8.42, 13.14), (79) 

(rc in kpc, rsoo hr Mpc, Mgas 500 hr IO^^Mq, M* 5 oo in IO^^Mq). The dif¬ 
ferent compoirents of the clnster halo iirtegrate to masses (np to distance r) 
docnmeirted in fig. 

The scalar field halo (SF) of the total baryonic mass psfi contribntes the 
lioir share to the total transpareirt energy. The SF of the galaxy system psf 2 
carries abont as mnch energy as the net phantom compoireirt Pphi in the 
baryceirtric rest system of the clnster. It snrpasses the gas mass close to the 
reference radins rsoo- Fignre shows the fast increase of the gravitational 
mass of the scalar field halo of the galactic system. 

If we add all baryoiric and halo contribntions, the pictnre given in fig. 
[^emerges. It shows an enconragiirg agreemeirt of the observed valnes M 500 
nt rsoo = 1278 kpc with the prediction of the Weyl geometric halo model 
Mtoi(rsoo) (in the range of the observational errors and of model errors). 

+0 97 +0 79 

Mtot(r5oo) = 5.66-0(68 , Msoo = 6.55-o'-rs x 10^^ Mq (80) 
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Coma M(r) in (1 Mo) 



— M_{t} 

- M_{sf} 

- M_{ph1} 

- M_{sf2} 


Figure 2: Halo components of Coma cluster: transparent matter halo Mt = 
Mgf + Mp/ji, total scalar field (SF) halo Mgf, halo of freely falling galaxies 
Msf 2 , net phantom energy Mphi (in harycentric rest system). 

Empirical data (violet dot and bar): baryonic mass + Mgas 500 (dot) 

and M 500 (with error intervals) at rsoo = 1280 kpc. 


Coma M(r) in (1 Mq) 



r in kpc 


— M_{bar} 

— M_{gas} 

— M_{sf2} 
M_{star} 


Figure 3: Comparison of the contribution of the scalar field halo of the galax¬ 
ies Msf 2 with the baryonic mass for the Coma cluster (empirical data for 
Mbar 500 violct dot). 
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Coma M(r) in (1 Mq) 



- M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


Figure 4: Contribution of the baryonic mass Mhar, of the scalar field and the 
phantom energies Mgf, Mpfii to the transparent mass Mt and to the total 
mass Mtot = Mt + Mhar of the Coma cluster in the WST model. Model 
errors indicated at rsoo, ?’200 (black). Empirical data for Mi,ar (violet dot) 
and for the empirically determined total mass M 500 with error bars (violet) 
from ]28f . Additional empirical data at r 2 oo (yellowish) from )18^ . 


Model error bars have been estimated by varying the input data (79) in their 
respective error intervals 
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At the Abell radius r^ = 2.14 Mpc (in the convention of [18]) we get the 
model values 


Mtotir a) - n.9, Miens{rA)-7.7 xlO^^M© (81) 


for the dynamical total mass and the lensing mass (76), ^ Due to the net 


phantom energy the dynamical mass is about 20 % higher than the lensing 
mass. This is an effect by which the model can, in principle, be tested 
empirically and discriminated from others,^® The baryonic masses in the ft 


model are Mgas{rA) ~ 11.73x 10^^ Mq^ M^,{rA) ~ 1.83x 10^^ Mq. Thus not 
only the total mass Mtot{v 50 o) given by our model agrees with the empirical 


^®The reconstruction of star mass from observational raw data is a particularly delicate 
point. Depending on the assumptions on the stellar dynamics and the resulting data 
evaluation model “one can obtain up to a factor of 2 fewer stars” |28l p. 6]. For obtaining 
our model errors we allowed a variation in stellar masses by factors 0.5 and 2. 

^®The total (dynamical) transparent matter contribution is MtfrA) = 9.30, of which 
Msf 2 (rA) ~ 2.01 are due to the scalar field halo of the galaxy system. The net phantom 
energy amounts to MphiiXA) ~ 1-82 (all values in units of 10 ^'^Mq). 

®°The lensing data of mi seem consistent with such an observation, Mie„s{ri) = 
6.1 -3-5 X 10^^ Mq at ri = 2.5 Mpc, compared with Ma ~ 12.9 x 10^"^ Mq at 
va ~ 2.14 Mpc in nHj. But because of their exorbitant large error intervals these data are 
far from significant for our question. 
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value M 500 inside the error bounds; also and M 200 are reasonably well 
recovered. 

This is the case without assuming any component of particle dark matter 
besides the (real) energy of the scalar held and the (phantom) energy ascribed 
to the additional acceleration induced by the Weylian scale connection. 
This might still be a coincidence. In order to learn more about the question 
whether the hndings at the Coma cluster are exemplary or not, we have to 
consider the data of all 19 galaxy clusters, respectively the 17 of the reliable 
sub-ensemble. 


3.6 Halos and total mass for 17(+2) clusters of galaxies 


The mass values of the halo models for the 19 clusters are calculated as 


described in section 3.3 with the choice of fadeout functions beyond r = r2oo 
(see appendix). The results are documented in tables 2, 3 and hgs. (for 
Coma see hg. 4). The baryonic masse at r 2 oo {Mgas{r 20 o), of table 

3 result from an extrapolation of the empirical values at rsoo given in |28] to 
?’200 using the /3-models of section 3.4, item (1), for gas and star mass. 

For 15 clusters our model reproduces the total mass at rsoo correctly, i.e. 
inside the error margins of data and model. This is achieved without any 
further adjustable parameter, only on the basis of the parameters for the 
/3-model for baryonic mass and (rsoo), Tf*(rsoo)- The fading out func¬ 

tions do not intervene below r 2 oo- Moreover, for the majority of these, and 
paradoxically for all other four, the less precisely determined data at r2oo 
have overlapping Icr error intervals. This indicates a surprising agreement 
between the (theoretically derived) transparent halo and the empirically de¬ 
termined dark halo, Mt ~ Mdm- 

For 5 clusters, A85, A2255, A2589 and the outliers A2029, A2065, the 
error intervals of empirical data and model data do not overlap. For the first 
three of them (A85, A2255, A2589) the model predictions are consistent with 
the empirical data within doubled error intervals (2a range). Only the two 
outliers (A2029, A2065) he farther apartOtherwise the model data are 
in good agreement with an assumption of normally distributed statistical 
errors and with the assumption that the evaluation bias due to the use of 


the NFW prohle for dark matter (item (3) in section 3.2) does not shift the 
mass estimates outside the error intervals. 

All in all, the assessment of the WST-3L halo model has surprisingly 


well passed, in spite of the main caveat of item (3), section 3.2 The outcome 
found for Coma seems to be typical also for the other clusters. Moreover, 
the good agreement of the model with the data of the reliable sub-ensemble 
of 17 clusters supports the assumption stated in the last phrase (3), section 
3.2 (no large systematic errors due to data transfer from Einstein/ACHM 


®^A2029 has the surprising property that the empirical values for the total mass at rsoo 
surpasses the one at r 2 oo, M 500 > M 2 oo. 
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Table 2. Empirical values (M 5 oo,M 2 oo) and model values (Mtot) for 
total mass at rsoo; ^200 


Cluster 

^500 

Mtot (r 500 ) 



^"200 

0 

to 

0 

0 

-17200 

Coma 

1278. 

5.32 

+ 0.83 

- 0.58 

6.55 

+ 0.79 

- 0.79 

2300. 

13.19 

+ 2.36 

- 1.65 

13.84 

+ 1.49 
- 1.41 

A85 

1216. 

4.67 

+ 0.51 

- 0.36 

6.37 

+ 1.00 
- 1.00 

1900. 

9.70 

+ 1.09 
- 0.77 

7.71 

+ 0.8 

- 0.74 

A400 

712. 

1.26 

+ 0.25 

- 0.16 

1.83 

+ 0.39 

- 0.39 

1093. 

2.62 

+ 0.55 

- 0.35 

1.48 

+ 0.21 

- 0.18 

IIIZw54 

731. 

1.33 

+ 0.34 

- 0.27 

1.91 

+ 0.58 

- 0.58 

1350. 

3.09 

+ 1.30 
- 1.02 

2.81 

+ 2.74 

- 1.10 

A1367 

893. 

1.82 

+ 0.29 

- 0.19 

1.76 

+ 0.27 

- 0.54 

1529. 

4.28 

+ 0.81 

- 0.62 

4.06 

+ 0.45 

- 0.40 

MKW4 

580. 

0.54 

+ 0.09 

- 0.06 

0.50 

+ 0.14 

- 0.14 

857. 

1.10 

+ 0.19 

- 0.13 

0.71 

+ 0.07 

- 0.06 

ZwC1215 

1098. 

3.74 

+ 0.47 

- 0.32 

4.93 

+ 0.98 

- 0.98 

2093. 

9.32 

+ 1.49 
- 1.07 

10.37 

+ 3.51 

- 2.62 

A1650 

1087. 

3.42 

+ 0.63 

- 0.50 

3.44 

+ 0.66 

- 0.66 

2150. 

9.40 

+ 2.87 

- 2.31 

11.14 

+ 5.77 

- 3.46 

A1795 

1118. 

3.40 

+ 0.40 

- 0.26 

3.41 

+ 0.63 

- 0.63 

2136. 

9.32 

+ 1.11 
- 0.74 

10.99 

+ 2.26 

- 2.09 

MKW8 

715. 

0.86 

+ 0.18 

- 0.14 

0.62 

+ 0.12 

- 0.12 

1279. 

2.35 

+ 0.67 

- 0.54 

2.38 

+ 1.04 
- 0.59 

A2029 

1275. 

6.39 

+ 0.61 

- 0.43 

14.7 

+ 2.61 

- 2.61 

2286. 

15.83 

+ 1.57 
- 1.12 

13.42 

+ 2.43 

- 2.26 

A2052 

875. 

1.63 

+ 0.26 

- 0.17 

1.39 

+ 0.28 

- 0.28 

1250. 

2.94 

+ 0.47 

- 0.32 

2.21 

+ 0.06 

- 0.08 

MKW3S 

905. 

1.90 

+ 0.27 

- 0.18 

1.45 

+ 0.34 

- 0.34 

1450 

3.81 

+ 0.60 

- 0.40 

3.46 

+ 0.36 

- 0.34 

A2065 

1008. 

3.92 

+ 0.72 

- 0.66 

11.18 

+ 1.78 
- 1.78 

2450. 

11.90 

+ 8.20 

- 4.61 

16.69 

+ 21.34 

- 6.73 

A2142 

1449. 

7.15 

+ 0.68 

- 0.50 

7.36 

+ 1.25 
- 1.25 

2364. 

15.29 

+ 1.53 
- 1.14 

15.03 

+ 3.9 

- 2.64 

A2147 

1064. 

3.32 

+ 0.53 

- 0.40 

4.44 

+ 0.67 

- 0.67 

1450 

5.81 

+ 1.14 
- 0.90 

3.46 

+ 1.17 
- 0.74 

A2199 

957. 

2.26 

+ 0.38 

- 0.28 

2.69 

+ 0.42 

- 0.42 

1621. 

4.98 

+ 0.93 

- 0.67 

4.81 

+ 0.37 

- 0.36 

A2255 

1072. 

3.94 

+ 0.44 

- 0.31 

7.13 

+ 1.38 
- 1.37 

2271. 

12.27 

+ 1.87 
- 1.37 

13.32 

+ 1.44 
- 1.19 

A2589 

848. 

1.92 

+ 0.33 

- 0.22 

3.03 

+ 0.75 

- 0.75 

1471. 

4.60 

+ 0.84 

- 0.58 

3.58 

+ 3.86 

- 1.54 


Model values MtotirNoo) empirical values Mvoo 10^^ Mq, 
fNOO (empirical) in kpc (N = 1, 2) 


to the WST-3L framework). But still we cannot exclude the possibility of 
cancelling between model errors and data transfer errors. Thus we have 
only found empirical support for the conjecture that, on the level of galaxy 
clusters, the observed dark matter effects encoded in M^m may solely be due 
to the combined impact of the halo Mgf of the scalar field and of the scale 
connection, Mphi (cf. @): 

Mdm ~ Mgf + Mph I = Mt (82) 


35 









Table 3. Model values for halo and baryonic masses at r 2 oo 


Cluster 

Mt 

Msf 

Msf2 

Mphi 

Mgas 

M* 

f* 

ft 

Coma 

11.15 

6.44 

1.73 

4.71 

1.77 

0.276 

0.16 

6.3 

A85 

8.03 

4.52 

1.01 

3.51 

1.54 

0.140 

0.09 

5.2 

A400 

2.27 

1.36 

0.45 

0.91 

0.26 

0.084 

0.32 

8.7 

IIIZw54 

2.77 

1.65 

0.54 

1.11 

0.25 

0.080 

0.32 

10.9 

A1367 

3.77 

2.21 

0.65 

1.56 

0.42 

0.089 

0.21 

8.9 

MKW4 

0.99 

0.58 

0.18 

0.40 

0.09 

0.022 

0.25 

10.9 

ZwC1215 

8.00 

4.55 

1.11 

3.45 

1.19 

0.137 

0.12 

6.7 

A1650 

8.15 

4.69 

1.24 

3.45 

1.10 

0.162 

0.15 

7.4 

A1795 

8.04 

4.59 

1.14 

3.45 

1.15 

0.14 

0.12 

7.0 

MKW8 

2.12 

1.24 

0.36 

0.88 

0.20 

0.039 

0.20 

10.8 

A2029 

12.83 

7.15 

1.47 

5.68 

2.83 

0.203 

0.07 

4.5 

A2052 

2.58 

1.50 

0.43 

1.07 

0.31 

0.059 

0.19 

8.3 

MKW3S 

3.35 

1.95 

0.55 

1.40 

0.39 

0.072 

0.18 

8.5 

A2065 

10.28 

5.80 

1.32 

4.48 

1.48 

0.142 

0.10 

6.9 

A2142 

12.55 

6.95 

1.35 

5.60 

2.59 

0.159 

0.06 

4.8 

A2147 

4.83 

2.77 

0.71 

2.06 

0.87 

0.118 

0.14 

5.5 

A2199 

4.36 

2.52 

0.68 

1.84 

0.55 

0.088 

0.16 

8.0 

A2255 

10.35 

5.84 

1.33 

4.51 

1.76 

0.166 

0.09 

5.9 

A2589 

3.99 

2.33 

0.68 

1.65 

0.52 

0.105 

0.20 

7.7 


Mass values in 10^^ Mq, /* = ft = ]^(r 2 oo), 

r 2 oo see tab. 2 
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A85 M(r) in (1 Mq) 






- M_{tot} 

— M_{t} 

— M_{sf} 
- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 
- M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


A400 M(r) in (1 Mq) 






- M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


Figure 5: Halo models for clusters 2 - 9 in tab. 1.: total mass Mtot (black 
line) with model error bars at rsoo; ^200; transparent matter halo Mt con¬ 
stituted by scalar field halo Mgf 2 net phantom halo (in barycentric rest 
system of cluster) Mphi and baryonic mass (gas and stars) M^ar- Empir¬ 
ical data for the total mass with error intervals at rsoo (violet) from ^281/ . 
Additional empirical data at r2oo (yellow) from fl8^ . 
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MKW8 M(r) in (1 Mq) 







- M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


A2029 M(r) in (1 Mq) 





- M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 



— M_{tot} 

— M_{t} 

— M_{sf} 

— M_{ph1} 

— M_{bar} 


— M_{tot} 

— M_{t} 

- M_{sf} 

- M_{ph1} 

— M_{bar} 


Figure 6: Halo models for clusters 10 - 19 in tab. 1. Description see fig. 
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3.7 Comparison with TeVeS and NFW halos 


It is surprising that in the WST-3L approach the total amount of observed 
dark matter Mdm seems to be explained by the energy of the scalar field 
and the phantom halo, Mt = Mgf + Mp^i, M^m ~ Mt- No missing mass is 
left. In usual relativistic MOND approaches this is not the case for clusters, 
although it is essentially so for galaxies [9]. Based on a study of about 
40 galaxy clusters, R Sanders has proposed the hypothesis of a neutrino 
component “between a few times 10^^ and 10^^ \ mostly concentrated 

close to the center of the cluster, in a region up to twice the core radius, 
supplementing the baryonic mass and the phantom energy of the TeVeS 
model m p. 902], In our approach, this hypothesis is unnecessary. Where 
does this difference arise from? 

A model calculation for the Coma cluster, evaluating (57) for the transi- 


2 (53) and the corresponding t '2 (54) used by 

1.8 X 


tion function fj .2 = x{l + x‘^)~ 

Sanders in |2n) . shows that the neutrino core had to be tuned to 
10^^ Mq in order to give agreement with the empirical value M^m(^50o) = 
-^500 ~ Mi)ar{r^QQ) Ri 4.7 X 10^^ Mq, where in this framework “rsoo”, “’’ 200 ” 
are to be defined by a formal convention with regard to a (fictitious) NFW 
halo. 

Table 3 shows the amount of scalar field energy up to a radius r ~ r 2 oo 
in the WST model. It varies between about 2 x 10^^ and 9 x 10^^ Mq, i.e., 
roughly in the range found necessary by R. Sanders for the (hypothetical) 
neutrino halo. Moreover, a comparison of the transition functions Hw{t) ( |55[ ) 
and lJi 2 {x) shows that the gross phantom energy pt of the WST approach is 
larger than in an /i 2 -MOND model |23| . Roughly half of the missing mass 
of Sanders’ model is covered by this effect, the other half is due to the scalar 
field halo of the system of galaxies up to r2oo- 

A comparison of the two MOND-like approaches is given in fig. Here 
one has to keep in mind that the TeVeS-/r 2 model has a free adaptable pa¬ 
rameter (mass of the neutrino core), while the WST has not. The general 
profiles of the “dark matter” halos of both models are similar. The TeVeS-/i 2 
mass starts from a higher socle because of its neutrino core; the WST trans¬ 
parent mass starts from a lower initial value but rises faster because of the 
increasing contribution of the scalar field halo of the galaxies. 

On the other hand, the best known profile for dark matter distribution, 
used in most structure formation simulations, is the NFW halo (Navarro/ 
Frank/White). Its profile is 


= A n + jL)2 ’ 

with density parameter po and core radius Tc (at which the density has 
reduced to half the reference value). For a first comparison of the interior 
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M_{dm} in (1 Mq) 



- M_{tWST} 

— M_{dm NFW} 

— M_{phm TeVeS} 


Figure 7: Comparison of dark/transparent/phantom mass halos for Coma 
in NFW, WST and TeVeS models, free parameters of halos for NFW and 
TeVeS {^2 with neutrino core) adapted to mass data (black error bars) at 
^500 = 1280 kpc 


halos we take Tc ~ 180 h~^ kpc {h = 0.7), following |13j. ^ and determine 
the central density parameter po such that the total integrated mass Mpfpw 
assumes the empirical value of m ^500 = 1289 A:pc (fig. [^. The error 
intervals for rsoo give upper and lower model values for po and corresponding 
model error bars (red) for Mi^Fw{'f' 20 o)- This version of the NFW halo 
satisfies one set of empirical data by construction (at rsoo); our main interest 
thus goes to the other empirical data set available at r 2 oo = 2300 kpc. The 
NFW model error interval overlaps with the empirical error bar at r2oo and 
with the error interval of the WST model. 

There is a conspicuous difference between the mass profiles of the NFW 
halo on one side and the halo profiles of WST or /r 2 -TeVeS on the other. 
More and precise empirical data of the interior halos of galaxy clusters ought 
to be able to discriminate between the two model classes. An empirical 
discrimination between the two MOND-like approaches would need more, 
and more precise, profile data. At the moment the Weyl MONDlike model 
survives the comparison fairly well, even though it has no free parameter 
which would allow to adapt it to the halo data. 


3.8 A side-glance at the bullet cluster 


At the end of this section let us shed a side-glance at the bullet cluster 
1E0657-56. It is often claimed that the latter provides direct evidence in 
favour of particle dark matter and rules out alternative gravity approaches. 
Our considerations show that this argument is not compelling. The energy 
content of the scalar field halos of the colliding clusters endows them with 

studies of the exterior halo of the Coma cluster, Geller e.a. have found fit values 

-b0.03 -1-0.029 -b0.035 , 

0.182 -0 03^0.167-0 029^0.192-0 035 in units h ^Mvc for |13|. 
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inertia of their own. The shock of the colliding gas exerts dynamical forces 
on the gas masses only, not directly on the scalar field halos. During the 
encounter the halos will roughly follow the inertial trajectories of their re¬ 
spective clusters before collision, and they will continue to do so for a while. 
It will take time before a re-adaptation of the mass systems and the respec¬ 
tive scalar field halos has taken place. Clearly the MOND-approximation 
is unable to cover such violent dynamical processes. It describes only the 
relatively stable states before collision and ~ in some distant future - af¬ 
ter collision. But a separation of halos and gas masses for a (cosmically 
“short”) period is to be expected, just like in the case of a particle halo with 
appropriate clustering properties. 

For the time being, the cluster 1E0657-56 does not help to decide between 
the overarching alternative research strategies, particle dark matter or alter¬ 
native gravity. It may be able to do so, once the dynamics of gas and of the 
halos has been modeled with sufficient precision in both approaches. Then a 
proper comparison can be made. But that is an overtly complicated task. It 
seems more likely that other types of empirical evidence will offer a simpler 
path to a differential evaluation of the two strategies and help clarifying the 
alternative. 


Discussion 


We have analyzed a three-component halo model for clusters of galaxies, 
consisting of 

(i) the scalar field energy induced by the overall baryonic matter in the 
barycentric rest system of a cluster (under abstraction of the discrete 


structure of the star matter clustering in freely falling galaxies) (69), 
with integrated mass equivalent Mgfi, 

[ii) an additional contribution to the scalar field energy, forming around 


the freely falling galaxies (71), integrating to Msf 2 


{in) 


and the phantom energy of the total baryonic mass in the barycentric 
rest system of the cluster, due to the additional acceleration of the 


Weylian scale connection in Einstein gauge (70) with mass equivalent 

Mphi- 

The first two components add up to a real energy content of the scalar field 
with mass equivalent Mgj = Mgfi + Msf 2 , the third one, Mphi, arises from 
a theoretical attribution in a Newtonian perspective and has fictitious char¬ 
acter. The mass equivalent of the integrated energy components combines 


to a total dark-matter-like quantity Mt = Mgf i -|- Mgf 2 -|- Mph 1 (73). 

All three components arise from gravitational effects of the cluster’s bary¬ 
onic mass Mbar = Mgas + Mstar in the framework of a Weyl geometric scalar 


41 






tensor theory of gravity 0 , with its scale connection as the specihc difference 
to Riemannian geometry 0, (|4) A second speciality of the theoretical 
framework is the cubic kinetic term of the Lagrangian (28), analogue to the 
AQUAL approach but in scale covariant form. Observable quantities are 
directly given by the model in Einstein gauge 0. 

The total dynamical mass of the model, Mtot = Mi)ar + Mt (75) has been 
heuristically confronted with the empirical values for of 17+2 galaxy clusters 
given in 125], |27], |28| . complemented by data from |18) . The problem of 
data transfer between different theoretical frameworks (in particular between 
Einstein gravity /ACDM and WST-3L) leads to a certain caveat with regard 
to an uncorrected taking over of the values for the total mass. But it does not 
seem to obstruct a meaningful hrst comparison of the WST halo model with 
available mass data of clusters collected in the Einstein/ACHM framework. 

The result of this comparison shows a surprisingly good agreement of 
the total mass predicted by the model Mtot, on the basis of data for the 
baryonic mass components, with the observed total mass M 500 (measured at 
the main reference distance ^ 500 )- Moreover the model shows an acceptably 
good agreement with additional empirical values at the distance r 2 oo given in 
|18) (determined on a slightly less rehned data basis and evaluation method). 
For 14 clusters the model predicts values for Mtot{fbOo) with error intervals 
(due to the observational errors for the baryonic data) which overlap with the 
empirical error intervals of M^qq. The Coma cluster is among them. Three 
clusters have overlapping error intervals in the 2a range. The remaining two 
are outliers and have been identihed as such already in EZj. 

In the result we have found empirical support for conjecturing that the 
observed dark matter at galaxy cluster level may he due to the transparent 
halo of the scalar field and the phantom halo of the scale connection of WST 

Details for the constitution of the total transparent matter halo from 
its specihc components (i), (u), (Hi) have been investigated for the Coma 
cluster (section |3.5[ ) . They seem to be exemplary for the whole collection of 
galaxy clusters. A particular feature of the model is the sealar field energy 
formed in the interspaces between the galaxies. Its integrated energy contri¬ 
bution surpasses the baryonic mass between 1 and 1.5 Mpc (see hg. [^and 
table 3, col. 4). It is crucial for this model’s capacity to explain the total 
dynamic mass on purely gravitational grounds, without any additional dark 
matter component. An overall comparison with R. Sanders’ TeVeS-MOND 


model for galaxy clusters and the NEW halo is given in section 3.7 


At the moment the Weyl geometric scalar tensor model with a cubic term 
in the kinetic Lagrangian of the scalar held fares well in all the mentioned re- 


Weyl geometric scalar tensor theory the scalar field (j) is the new dynamical variable, 
while the scale curvature dtp = / vanishes. The latter would be an even more striking, 
and even irritating, difference to Riemannian geometry; in integrable Weyl geometry it 
plays no role. 
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spects. It would be very helpful if astronomers decided to evaluate old or new 
raw data in the framework of WST. That could lead to an empirical discrim¬ 
ination of the different models. But already independent of the outcome of 
such a revision, the scalar field cj) and the scale connection of WST have a 
remarkable property from a theoretical point of view. They complement the 
classical Einstein-Riemannian expression for the gravitational structure, the 
metric field by a feature which carries a proper energy-momentum ten¬ 
sor (35), (36). The energy momentum of the scalar field plays a crucial role 
for the constitution of the transparent matter halo. In the present approach 
it seems to express the self-energy of the extended gravitational structure. It 
remains to be seen whether this is more than a model artefact. 


5 Appendix 

5.1 Energy component of the scalar field 

In scalar field (Einstein) gauge = — fiifu = 

and Dyffi = 2(j)Di^(j) = 2((P‘duU}. Simliarly 


Moreover, 


D^d^'u) = V^d'"oj-2ip^d''uj =Pfj,d''oj+5''^ip\d uj-g^xtp^d oj+2d^u]S'"uj 

leads to 

Dfj_D^(j)^ = 2cj)^{2dfj_ujd’^uj -\-2dfj_ujd’^uj — d’j^dxiod^uj) 

= 2cj)^{fi7fj_d^uj->r 4:d^ujd’^uj — S’j^dxLvd^io). 

Thus 


DxD^ct? = 2fi^ gVxd^uj 


(84) 


Following (35), (36), (33), the energy momentum of the scalar field, 

0 = (r(^) + , 


becomes 


= -(SttG)-^ (^2gVxd^uj -Sdxojd^eo - ‘^do\Vu\^ + g 


= eDgD,fi^-2 


dL^ 

dgU^ 


(SttG) ^ + (1 - do\Vuj\)dij,ujdgUJ - dx^xid^u g^y^ . 
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In the static case the hrst two terms of T, 
component of the scalar held turns into 


(ll) 


fJ,U 


Ti = - dxujd^uj- 


vanish and the total energy 

+ . (85) 


Assuming conditions under which the terms of cosmological orders of mag¬ 
nitude and those of order iVwp can be neglected, we arrive at 




like in equation (46) of the main text. In the central symmetric case ( |44[ ) 
implies uj{r) = y/GMao Inr. In spherical coordinates Vw = (0, ^ Q) 

and 


72 V7 oA . • VGMdo 


=„Vxd^ui = 


( 86 ) 


Because of diudiu = 


__— ^ ^ the approximation (46) is justihed. 


5.2 Remarks on the numerical implementation 

The calculations in sections 3.2ff. have been implemented in Mathematica 
10 and run on a PC. Integrations of the mass values have been realized by 
numerical interpolation routines in distance intervals of 100 kpc. A compar¬ 
ison with rehned distance intervals 10 kpc showed differences at the order of 
magnitude 10“^ of the respective values, thus below the rounding precision. 

The fading out for the scalar held and phantom halos beyond r2oo has 
been modeled by the cubic expression: 


f{x, A, B) = x{x; -oo. A) + 


1 


1 + 


x—A 

B 


3 

x(x;A,oo ), 


(87) 


with xix',CL,b) the characteristic function of the interval [a,b]. The fading 
out of /(x. A, B) starts at A and declines to ^ A + B. In our case we start 
the fading out close to the virial radius, A = l.lr2oo ^-nd set B = 0.5r2oo- 
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